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Simulation studies are performed to assess the sensitivity of a model- 
independent analysis of the flavour-tagged decays D° — > K^tt + tt" and 
D° — > K§K + K~ to mixing and CP violation. The analysis takes as input 
measurements of the D decay strong-phase parameters that are accessible 
in quantum-correlated D D pairs produced at the -?/>(3770) resonance. It 
is shown that the model-independent approach is well suited to the very 
large data sets expected at an upgraded LHCb experiment, or future high 
luminosity e + e~ facility, and that with 100M K^tt + it~ decays a statistical 
precision of around 0.01 and 0.7° is achievable on the CP violation pa- 
rameters r C p and acp, respectively. Even with this very large sample the 
systematic uncertainties associated with the strong-phase parameters will 
not be limiting, assuming that full use is made of the available -^(3770) 
data sets of CLEO-c and BES-III. Furthermore, it is demonstrated that 
large flavour-tagged samples can themselves be exploited to provide infor- 
mation on the strong-phase parameters, a feature that will be beneficial 
in the measurement of the CKM angle 7/^3 with B~ — > DK~ decays. 



1 Introduction 



In the last few years measurements have been performed at the 5-factories and the 
Tevatron which, when taken together, reveal the presence of mixing in the D°-D° 
system [l}j6]. Improved sensitivity to this phenomenon will come from the LHCb ex- 
periment, which has already published a first study with early data [7]. Interest in the 
charm sector has been raised still further with the recent announcement of evidence 
for CP violation in time-integrated D° — > K + K~ and D° — > tt + tt~ ~ decay^Hg. This 
effect, if confirmed by future measurements, is almost certainly attributable to direct 



CP violation 10 . A natural next step in the experimental programme is to intensify 
the search for indirect CP violation, which manifests itself through mixing-related 
observables. In the Standard Model (SM) such CP violation is expected to be very 
small and so any non-zero effect would be indicative of possible New Physics (NP) 



contributions 11 



A particularly powerful decay mode for flavour-tagged time-dependent D°-D° 
studies is the channel D° — > K^tt + tc^, where the interference between intermediate 
resonances provides sensitivity to both the magnitude and sign of the mixing pa- 
rameters, and allows for probes of indirect CP violation. The analysis of this mode 
was pioneered by CLEO [12], and higher precision results have been published by 
Belle (4] and BaBar |5|, based on samples of around 0.5 M and 0.7M events respec- 
tively. Studies based on larger data sets can soon be expected from LHCb. In time, 
sample sizes of O(10 8 ) events may become available from future experiments, such as 



Belle-II 13 , SuperB [14] and the LHCb upgrade 15 16 . It is important, therefore, 
to understand what sensitivity to CP violation can be attained with large data sets of 
D° — > K^tt + tt~ decays, and what challenges such a measurement will be confronted 
with. 

The D° — > fCg7T + 7r _ analysis is performed in the Dalitz space of the three-body 
final state. Measurement of the mixing and CP violation parameters requires a 
good understanding of the decay structure within this Dalitz space. The established 
method is to develop an amplitude model of the D° decay on flavour-tagged data, and 
then exploit this model in an unbinned time-dependent likelihood fit. The drawback 
to this approach is that the results have a systematic uncertainty that is associated 
with the assumptions of the model. The size of the uncertainty is difficult to assess, 
and the degree to which it can be decreased with larger data samples is unknown. An 



alternative strategy, first proposed in 17 , is to perform a binned, model-independent 
analysis, where the necessary inputs on the D° decay come from auxiliary measure- 
ments performed on quantum-correlated D D decays at the ^(3770) resonance. In 
this case the systematic uncertainties associated with the model are replaced with 
measurement uncertainties coming from these external inputs. No data analysis has 
yet been presented which follows this model-independent approach. 



^Unless stated otherwise, the charge-conjugate mode is implicit throughout this article. 
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Identical considerations apply to the channel D° — > K§K + K~ . Although sup- 
pressed by a factor of around six compared to D° — > K^-k + -k^, this decay can be 
exploited in an identical manner for mixing and CP violation studies. A model- 
dependent analysis has been performed by BaBar [H], and Belle have studied the 
time-evolution of the CP-odd component of the decay [6). The model-independent 
study of the full Dalitz space is also an appealing analysis option for this decay mode, 
but one that has not yet been pursued. 

This paper presents the results of simulation studies performed to investigate the 
sensitivity of the model-independent approach to the decays D° — > fTg7r + 7r _ and 
D° — > K§K + K~ . The statistical precision of a binned analysis is assessed and com- 
pared with that of a model-dependent, unbinned study. The sensitivity is calculated 
both for the sample sizes that currently exist, and for those which will become avail- 
able over the coming years. When considering future large data sets, the point at 
which the analysis will become limited by the existing measurements of the D meson 



decay properties from the ?/>(3770), available from the CLEO-c experiment 18 , is de- 
termined, together with what improvement in the knowledge of these properties will 
be required. Various fit strategies are explored which make different use of external 
inputs. It is shown that other useful results can be obtained from the analysis, aside 
from the mixing and CP violation parameters. 

Model-independent studies of mixing and CP violation with other D° decay modes 
are possible (e.g. see Ref. (19)), all of which can benefit from measurements of strong- 



phase parameters and related properties in -^(3770) data 20 



2 Formalism 

In the neutral D meson system the mass eigenstates, -D 12 , are related to the flavour 
eigenstates D° and D° as follows: 

\D 1>2 ) = p\D°) ± q\D°) (1) 

where the coefficients satisfy \p\ 2 + \q\ 2 = 1 and 

r CP e iacp = -. (2) 
V 

Indirect CP violation occurs if tqp ^ 1 and/or acp ^ 0. Charm mixing is conven- 
tionally parameterised by the quantities xn and yn, defined as 

m 2 - Mi r 2 - r\ 

%d = ^ , Vd = 2r (3) 

where Mi j2 and I\ 2 are the mass and width of the two neutral D meson mass eigen- 
states, and T the mean decay width of the mass eigenstates. A global analysis of all 
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yields the results x D = (0.631^)%, y D = (0.75 ±0.12)% 

+9.4\o 



lO.l^g g)°. The values of the latter two parameters are 



current measurements 10 
r CP = 0.89+°;^ and a CP = 
consistent with the null CP violation hypothesis. 

The time evolution of neutral D° mesons proceeds as follows: 

\D°(t))= g + (t)\D°) + ?- g _(t)\D ), 

V 

\D (t)} = Pg4t)\D°}+ g + (t)\D°), 



(4) 
(5) 



where 



-i(M 1 -iT 1 /2)t _|_ e -i(M 2 -ir 2 /2)t\ 



(6) 



Consideration is now given to the specific case in which the D° decays to the final 
state .Kg7r + 7r~ (the formalism presented is identical for K$K + K~). The final-state 
particles are labeled as follows: 1 for K®, 2 for 7r + and 3 for 7r~. The Dalitz plot 
decay amplitude density of D° — > K^-n + n^ is defined as 

A D o(m 2 12 ,m 2 13 ) = a 12il3 e^ 2 - 13 , (7) 

where m\ 2 an d specify coordinates on the D° — > Xg7r + 7r~ Dalitz plot, 012,13 is 
the modulus of the amplitude and #12, 13 is a CP-conserving strong-phase. The decay 
amplitude density is independent of the time evolution of the D° system. Assuming 
there is no direct CP violation in decays to the final state in question, the equivalent 
decay amplitude density for D° decays is 



ApAl in 



2 2 \ 
12' m 13) 



A D a(ml 3 ,ml 2 ) 



oi$i3, 12 



Ol3,12e 

The time-dependent probability density in which mixing and indirect CP violation 
are allowed is 

2 



V D o{m 2 l2l m 2 l3 ,t) 



g + (t)A D o(m 2 12 ,m 2 13 ) + -g^(t)A 50 (m 2 12 ,m 2 13 ) 
2 



T\g + (t)\ 2 a 2 



12,13 



+ r 



\g-(t)\ 2 a 



2„2 
13,12 



n* (+\ n n <5l2, 13+^13,12) 

~9+\ t )9-V') a 12,13 <2l3,12 e 



V 



g+{t)g*_{t) a 12A3 a 13 , 12 e ti5 ^- 5 ^ 



(9) 



which, ignoring terms of O^x 2 -)), 0(xdVd) and 0{y 2 D ), is equal to 



V D o(m 2 12 ,m 2 13 ,t) = Te 



-n 



l 12,13 



r C p 012,13 013,12 rt {y D cos(<5 



12,13 
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13,12 — Oicp) 



+ x D sin(5i2,i 3 - #13,12 - a C p)} 



(10) 
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Integrating the probability density between times t a and U yields 

V D o(m 2 12 ,m 2 13 ,t) dt = 

" (e~ rt « - e- rt »)a 2 12A3 + [T(e~ r H a - e~ rt n b ) + (e^ - e^)] 

x ^012,13^13,12 { Vd [cos(5i2,i3 - £13,12) cos(a CP ) + sin(5 12i i3 - #13,12) sxa(a C p)) 
+ x D [sin(<5i 2j i3 - £13,12) cos(a CP ) - cos(£i 2 ,i 3 - £13,12) sin(a CP )]}. (11) 

The formalism for the time evolution of the D° decay probability is identical, except 
rep is replaced by l/rcp and a>cp is replaced by —a>cp- 

If a model is available to describe the variation of amplitude across the Dalitz plot 
then Eqn. [IT] may be used in a maximum likelihood fit, unbinned in Dalitz space, to 
determine the mixing and CP violation parameters in the analysis of a given data 
sample. If however a model-independent approach is pursued then it is necessary 
to partition the Dalitz plane. A total of 2N regions are defined, symmetric under 
the exchange m\ 2 -H- m\ 3 and labelled from —N to +iV (excluding zero). Then the 
population in bin % is given by 

/ / V D o [m\ 2 ,m\ 3 ,t)dt dm 2 l2 dm 2 l3 = 

Ji Jt a 

n{{e- vt « - e- Ttb )Ti + [r(e" r ^ a - e' Tt H b ) + (e^ - e"^)] 

x [r C p\/T{rZ i {y D [ciCOs{acp) + s { sm(a CP )] + x D [si cos(a C p) - q sin(a C p)])}} 

(12) 

where n is a normalisation factor and Tj, q and Sj are defined as follows: 

Ti = Ja\ 213 dm\ 2 dm 2 13 , (13) 
y 012,13 «i3,i2 cos(£i 2; i3 - £13,12) dmf 2 dm\ 3 , (14) 



\[W- 

r—— / «i2,i3 013,12 sin(£i2,i3 - £13,12) dm? 2 dm\ 3 . 
\/TiT_i Ji 



(15) 



These expressions make use of the relations c_j = q and s_j = — s$. In the limit of 
no mixing or CP violation Tj is proportional to the population in bin % coming from 
decays of D mesons of a known flavour. The quantity q (sj) is the cosine (sine) of 
the average strong-phase difference between D° and D° decays within bin i, weighted 
by the decay rate. 

As external measurements are available for q and Sj (see Sect. [3]), and the values 
of Tj may also be obtained from external data sets, or indeed left as free parameters 
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in the fit (see Sect. 4.1), it is possible to use the observed population of the Dalitz 
and decay time bins to determine the mixing and CP parameters independent of any 
model assumptions. 



3 Available and future measurements of Cj and Si 

The CLEO collaboration have published results which exploit the quantum-correlated 
nature of ^(3770) — > DD decays to measure the Cj and coefficients for both 
D° K^iT + 7r~ and D° K^K + K~ (l8|. The study was motivated by the need 
to provide input for the model-independent measurement of the CKM unitarity an- 



gle 7/03 in B~ — > DK~ decays 21 22 . A data set of 0.8fb _1 is analysed and the 
parameters extracted from measuring the relative rates of various event classes, for 
example events containing both a neutral D decay to the signal mode and a neutral 
D decay to an even or odd CP eigenstate. Binning schemes are defined following the 
guidance of amplitude models developed at the B-factories with the goal of providing 
a partitioning which provides good statistical sensitivity in the 7/^3 analysis. The 
use of an amplitude model to define the binning introduces no model-dependent bias 
in the measurement of the q and Sj parameters. However if the model does not de- 
scribe nature well then the statistical power of the binning choice will be diluted with 
respect to expectation. 

In D° — > i^g7r + 7r~ decays the coefficients are measured in 8 pairs of Dalitz plane 
bins. Results are given for four choices of binning: two where the Dalitz space is par- 
titioned into equal intervals of the strong-phase difference, with the variation of the 



strong-phase difference assumed to follow the models developed by either BaBar 23 
('equal A8d BaBar') or Belle [24] ('equal A8u Belle'); and two where the binning is 
optimised to yield the best possible sensitivity to 7/^3 in a low background ('optimal 
binning') and higher background ('modified optimal binning') environment. For the 
equal interval strong-phase analysis the typical uncertainties on each q (s») measure- 
ment range from around 0.06 to 0.16 (0.11 to 0.21), depending on the bin. 

In D° K^K + K~ decays the q and Si coefficients are measured in an equal 



interval strong-phase binning, derived from the model presented in 25 . The lower 
statistics for this decay restrict the analysis to a fewer number of bins than is the case 
for D° — > K^7c +/ k~ . Results are presented for the Dalitz space partitioned into pairs 
of 2, 3 and 4 bins. With 4 pairs of bins the precision on q varies between 0.07 and 
0.36 for q, and 0.31 and 0.73 for Sj. 

The bin definitions for the equal A5d BaBar fTg7r + 7r~ scheme and that for K%K + K~ 
with 4 pairs of bins is shown in Fig. [Tj 

The knowledge of q and Si for both decays can be improved by future measure- 
ments. Already the BES III experiment has accumulated around 3 fb _1 at the ^(3770) 
resonance and a total integrated luminosity of 5 — 10 fb _1 is foreseen at this collision 
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Figure 1: Equal A5d BaBar binning for K^7t + tt (left) and for K$K + K (right). 



energy 26 . Future facilities, such as a r-charm factory in Novosibirsk 27 or the 
SuperB project operating at the charm threshold jT4] would provide substantially 
larger samples of up to 500 fb _1 . 

4 Fit studies 

In order to determine the projected uncertainties on the mixing and CP violation 
parameters in the model-independent approach a programme of 'toy Monte Carlo' 
studies is conducted. In each study events are generated in bins of Dalitz space and 



D° decay time according to Eqn. 12 Unless stated otherwise, the binning scheme 
adopted is the equal A5d BaBar binning with 8 pairs of bins for .K"g7r + 7r~ and the 4 
pairs of bins for K$K + K~ . For the results reported below xr> and yo and the CP 



violation parameters are set to their current world average values 10 . The number 
of events in each bin is smeared to account for statistical fluctuations. No background 
or detector effects are included. A binned maximum likelihood fit is then applied to 
determine the parameters of interest. Many samples are generated and fit, and in 
the case of successful convergence the fit parameters are recorded for each simulated 
experiment. The precision on each parameter of interest is determined by taking the 
standard deviation of the set of fit values. 

In general the conclusions of each study have a dependence on the sample size. 
Therefore three scenarios are considered: 

1. Data sets of a similar magnitude to those that have already been analysed at 
the 5-factory experiments. In the case of K^7t + ii^ this is 0.5M decays, and for 
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K^K+K- 0.1M decays. 



2. Data sets that may become available at LHCb after several years of operation. 
Sample sizes of 10M and 2M are assumed for J\g7r + 7r~ and K®K + K~ respec- 
tively. 

3. Data sets that may become available at an upgraded LHCb or a future high 
luminosity e + e~ experiment such as Belle-II or SuperB. Sample sizes of 100M 
and 20M are assumed for K"g7r + 7r~ and K$K + K~ respectively. 

In the case of 2 and 3 the numbers are indicative, and should not be interpreted as 
official estimates for any of the experiments mentioned. 



4.1 Fitting with external input for (q, Si) 

In the first set of studies the Cj and Sj parameters, as measured by CLEO-c, are 
taken as external input, and the mixing and CP parameters xd, Vd, rep and acp 
are determined from the fit. The results for 100M J\g7r + 7r _ decays are presented in 
Table [U 

In one set of experiments, denoted '(q, s») fixed', the values of q and Sj are fixed 
to be the CLEO-c central values in both the generation and fit; in another, denoted 
'(q, Si) smeared', the parameters used in the generation are smeared in a correlated 
manner according to their measurement uncertainties, but kept fixed to their central 
values in the fit. Comparison between these two scenarios allows the impact of the 
(q, S{) measurement uncertainties on the results for the mixing and CP violation 
parameters to be assessed. 

It is also possible to use external input for the Tj parameters. Results are available 



from CLEO-c 18 , but significantly more precise measurements can be made with the 
larger flavour tagged samples that are already available from other experiments. In 
the studies three scenarios are considered: fixing the values of Tj in the generation 
and the fit to those measured by CLEO-c ('Tj fixed'); smearing the values in the 
generation according to the measurement uncertainties that would be achievable with 
existing 5-factory data sets [4] , but keeping them fixed to their central values in the 
fit ('Tj smeared'); fixing the values of Tj in the generation but including them as free 
parameters in the fit ('Tj floating'). 

For sample sizes of a few million events and above it is found that the uncertainties 
on the current measurements of Tj limit the precision with which the mixing and CP 
parameters can be determined. An estimate of the systematic uncertainty on each 
parameter coming from this source can be obtained by subtracting in quadrature 
the results for 'Tj fixed' from those of 'Tj smeared'. This exercise yields errors of 
approximately 0.216(3)% for x D , 0.242(3)% for y D , 0.0483(7) for r C p and 3.65(5)° for 
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Table 1: Expected uncertainties on mixing and CP parameters for 100M K^-k^tt^ 
decays, in a variety of fit scenarios. The values in parentheses indicate the precision 
on each number. 



Scenario 


a(x D ) (%) 


*(lto) (%) 


a{r C p) 


cr(a C p) (°) 


(q, Si) fixed, Tj fixed 
(q, Si) fixed, Tj smeared 
(cj, Sj) fixed, T floating 
(q, Si) smeared, Tj floating 


0.0124(2) 
0.2168(29) 
0.0173(2) 
0.0757(10) 


0.0135(2) 
0.2424(32) 
0.0189(2) 
0.0886(11) 


0.0116(2) 
0.0497(7) 
0.0118(2) 
0.0262(3) 


0.740(10) 
3.721(49) 
0.744(10) 
1.188(15) 



olcp- Improved precision is obtained if the Tj values are left as free parameters in the 
fit. This approach is therefore considered as baseline for the subsequent discussion. 

The statistical uncertainties from the K^7i + 7i^ analysis, as estimated from the 
'(q, Sj) fixed, T floating' ensemble of experiments, for the three sample sizes under 
consideration are shown in Table [2] In addition are listed the estimated system- 
atic uncertainties in the binned analysis arising from the current knowledge of the 
(q, Sj) parameters. These errors are determined by subtracting in quadrature the 
mean uncertainties of the '(q, Sj) fixed, Tj fixed' event experiments from those of 
the corresponding '(q, Sj) smeared, Tj fixed' studies. Also shown are the statistical 
and model uncertainties of an unbinned model-dependent Belle analysis [I] conducted 
with around 0.5M decays. 

Table 2: Expected statistical and systematic uncertainties from knowledge of (q, Sj) in 
a binned K^n + 7T^ analysis for a variety of sample sizes. The values in parentheses in- 
dicate the precision on each number. Also shown are the statistical and model-related 
systematic uncertainties from an unbinned Belle analysis of 0.5M signal events (4l. 



Parameter 


Binned 


Unbinned, 0.5M 


Stat, 0.5M Stat, 10M Stat, 100M Syst (a, s<) 


Stat Syst model 


a(x D ) (%) 

<t(Vd) (%) 
<y{rcp) 
<7(a CP ) (°) 


0.251(3) 0.054(1) 0.017(0) 0.076(1) 
0.272(4) 0.061(1) 0.019(0) 0.087(1) 
0.175(2) 0.037(1) 0.012(0) 0.024(0) 
12.46(16) 2.42(3) 0.74(1) 0.88(2) 


30 +0 09 

u - ou -0.16 
U - Zd -0.08 

1°J° 0.08 

+16 +2 
-18 -4 



The results in Table |2] indicate that the statistical uncertainties in the binned 
analysis exhibit an approximate 1/y/N scaling, as expected. The 0.5M binned sta- 
tistical uncertainties are similar to those of the unbinned Belle analysis. Indeed, in 
some cases the binned uncertainties are slightly smaller than those from the unbinned 
study. It is assumed that this behaviour can be attributed to the idealised nature of 
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the toy Monte Carlo and the particular properties of the Belle data set. In principle 
the small loss of information that occurs in a binned study should lead to a slightly 
reduced statistical precision with respect to an unbinned analysis. Further discussion 
on this issue can be found in Sect. 14.41 

Comparison of the (q, Sj) systematic uncertainties with the statistical uncertain- 
ties associated with a data set of a given size allows the potential of the available 
CLEO-c measurements in the K^7i + ir~ mixing and CP violation analysis to be as- 
sessed. Of most interest are the CP violation parameters r C p and a>cp- Here the (c$, 

uncertainties are smaller than the statistical uncertainties for a 10M event data 
set. By interpolating the statistical precision it is found that the (q, s») uncertainty 
becomes dominant in data sets of around 25M and 75M events, for rep and acp 
respectively. Even in the case of 100M events it would be presumably be possible to 
make the (c*, Sj) uncertainty subdominant, merely by repeating the CLEO-c anal- 
ysis with the 3 fb _1 -^(3770) sample that BES-III now possesses. The current (q, 
Si) uncertainties are significantly smaller than the assigned model errors in the Belle 
analysis. Although larger data sets will allow for the K^tt + tc~ model to be refined, it 
is difficult to assess what improvements can be expected in the associated systematic 
uncertainties. 

The precision on the mixing parameters xp> and yp> more rapidly becomes limited 
by the current knowledge of (q, si). For 10M events these strong-phase uncertainties 
are dominant. They are of a similar size to the assigned model uncertainties in the 
unbinned Belle analysis. Analysis of the quantum-correlated D-D data set already 
available at BES-III would allow for the (q, Sj) uncertainty to remain subdominant 
for i^g7T + 7r~ sample sizes up to around 20M events. 

Table 3: Expected statistical and systematic uncertainties from knowledge of (q, s,) 
in a binned K®K + K~ analysis for a variety of sample sizes. The values in parentheses 
indicate the precision on each number. 



Parameter 


Stat, 0.1M 


Stat, 2M 


Stat, 20M 


Syst (q, Si) 


a(x D ) (%) 


0.807(11) 


0.168(2) 


0.066(1) 


0.330(4) 


<Vd) (%) 


0.546(7) 


0.102(1) 


0.033(0) 


0.178(2) 


o{.r CP ) 


0.266(4) 


0.073(1) 


0.029(0) 


0.068(1) 


cT(a CP ) (°) 


33.55(44) 


6.47(9) 


2.34(3) 


3.88(7) 



The statistical uncertainties that are expected for a binned analysis of K§K + K~ 
data, assuming the four-bin equal strong-phase binning, are shown in Table [3] As 
expected the precision is worse than in the iTg7r + 7r _ case. The uncertainties on the 
mixing parameters for the 0.1M scenario are similar to the values of 0.92% for xp> 
and 0.57% for y D that are obtained in an unbinned BaBar analysis [HI. The (q, Sj) 
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uncertainties from the CLEO-c analysis are dominant in the measurement of xp, and 
yr> for the 2M sample, but not so for r C p and acp- For 20M events the (q, si) 
uncertainties on r C p and acp are around twice the size of the expected statistical 
uncertainties. Exploitation of the already available BES-III -0(3770) data would allow 
for the (q, Si) uncertainties to be approximately halved. 



4.2 Fitting with external input for xd and yu 

An alternative strategy to that presented above is to consider r C p and a>cp as the 
principal parameters of interest, and to use external measurements to fix xd and yo 
in the fit. In this approach it is natural also to determine (cj, s«) from the fit, as in 
most scenarios this may be done with a precision that is similar to or better than that 
of the CLEO-c analysis. Improved knowledge of (q, Sj) is useful for the measurement 
of 7/03 using B~ — )■ D(K§h + h~)K~ decays. 

The free parameters are therefore rep, acp, (c?. ; Si), Tj, T and the overall number of 



events. The values of xd and yp> are fixed to their current world- average values 10 
in both the generation and the fit. The statistical uncertainties on rep and acp 
and the systematic uncertainties arising from the current knowledge of the mixing 
parameters, given separately for xp, and yp>, are shown in Table |4} The statistical 
precision is somewhat worse than that obtained following the procedure presented in 



Sect. 4.1 The systematic uncertainty from the current knowledge of (xp, Vd) does 
not dominate the measurement of rep until the sample size exceeds 10M decays, with 
the error associated with xp, being the most significant. For acp the measurement 
remains statistics-limited up to even larger data sets, making this an attractive and 



complementary alternative analysis strategy to that discussed in Sect. 4.1 



Table 4: Expected statistical and systematic uncertainties on the CP violation pa- 
rameters from knowledge of xp and yp> in a binned i^g7r + 7r~ analysis for a variety of 
sample sizes. The values in parentheses indicate the precision on each number. 



Parameter 


Stat, 0.5M 


Stat, 10M 


Stat, 100M 


Syst (x D ) 


Syst (y D ) 


c(Vcp) 


0.3129(44) 


0.0561(8) 


0.0170(2) 


0.0414(7) 


0.0234(5) 


v((*cp) (°) 


27.26(38) 


3.35(5) 


0.99(1) 


0.51(4) 


0.29(7) 



The uncertainties on q and Sj from a sample of 100M fCg7r + 7r~ events are pre- 
sented in Table [5} The results depend on how well the input parameters xd and yp 
are known, and so the study is repeated considering a variety of scenarios in which 
the uncertainty on both parameters is scaled from their current values down by a 
common factor. Also shown are the measurement uncertainties from the CLEO-c 
quantum-correlated analysis. It can be seen that even without any improvement in 
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the knowledge of the mixing parameters it is possible to fit q and Sj in certain bins 
with a precision that is similar or better than that of the CLEO-c measurements. 
If the knowledge of the mixing parameters improves by a factor of three or four 
in precision, for example through mixing analyses using alternative decay modes or 
through analysis of independent data sets, then it is possible to reduce significantly 
the uncertainties on the (q, Sj) parameters. 



Table 5: Fit uncertainties on and Sj for 100M K^n + ir~ decays when fixing (xjj, yo) 
and assuming a certain knowledge of these parameters with respect to the current 
precision. Also shown are the current measurement uncertainties on q and Sj from 
CLEO-c fl8l. 



Parameter 



(7(Ci) 

<x(c 2 ) 

cr(c 4 ) 
o-(c 5 ) 

cr(c 6 ) 

a(c 7 ) 
CT ( c s) 



cr(s 4 ) 

<r(s 7 ) 
o"(ss) 



Scale factor for current uncertainties on% and yr> 



1 



1/2 



1/3 



1/4 



l/oo 



0.125(2 
0.085(1 
0.068(1 
0.108(2 
0.159(2 
0.114(2 
0.065(1 
0.103(1 



0.102(1 
0.123(2 
0.225(3 
0.169(2 
0.055(1 
0.219(3 
0.236(3 
0.170(2 



0.087(1 
0.062(1 
0.064(1 
0.079(1 
0.091(1 
0.083(1 
0.056(1 
0.071(1 



0.102(1 
0.089(1 
0.154(2 
0.124(2 
0.052(1 
0.151(2 
0.161(2 
0.117(2 



0.079(1 
0.056(1 
0.064(1 
0.071(1 
0.071(1 
0.079(1 
0.056(1 
0.065(1 



0.100(1 
0.075(1 
0.118(2 
0.099(1 
0.051(1 
0.123(2 
0.119(2 
0.096(1 



0.076(1 
0.053(1 



0.062( 
0.0681 
0.0601 
0.0761 
0.0561 
0.062( 



0.098(1 
0.067(1 
0.099(1 
0.091(1 
0.051(1 
0.105(1 
0.094(1 
0.084(1 



0.072(1 
0.050(1 
0.064(1 
0.065(1 
0.048(1 
0.070(1 
0.056(1 
0.058(1 



0.097(1 
0.057(1 
0.076(1 
0.078(1 
0.049(1 
0.088(1 
0.065(1 
0.070(1 



CLEO-c 



0.055 
0.093 
0.161 
0.153 
0.062 
0.126 
0.163 
0.102 



0.107 
0.195 
0.138 
0.214 
0.134 
0.205 
0.154 
0.158 



As can be seen from Table [5] the statistical precision is better in some bins com- 
pared to others, and the measurement uncertainties from CLEO-c also vary from 
bin-to-bin. Therefore analysis of sample sizes of ~ 20 M decays would already be 
beneficial for improving knowledge of certain parameters. 

A similar analysis performed on K§K + K~ simulated data yields the same quali- 
tative conclusions. Again the precision on the CP violation parameters is inferior to 
that observed in Sect. 4.1, but it is possible to measure the (cj,Sj) parameters with a 
precision that is useful for future 7/^3 studies. 
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4.3 Fitting with no external inputs 



With a sufficiently large data set it is possible to determine the mixing, CP violation 
and (q, parameters in the fit. In practice this is achieved through a two-stage 



procedure. In the first stage the parameters of small magnitude in Equation 12 



namely xd, JJd, Ci and Sj, are fixed to their current world-average values and the 
other parameters are determined. In the second stage the parameters that were fitted 
in the first stage are fixed, and vice versa. 

When this study is performed for ensembles of 100M K^n + n~ decays uncertainties 
of 0.0644(10)%, 0.0790(12)%, 0.0244(4) and 1.762(27)° are obtained for x D , y D , r C p 
and acp respectively. These results are not competitive with those of the strategies 
that make use of external inputs. The precision for (q, si) in this study is on average 
around twice as poor as that obtained when constraining xe> and yu to lie within 
their current measurement uncertainties. 



4.4 Varying the number of bins 

It is interesting to consider the variation in the precision of the fitted parameters 
if the model-independent analysis is repeated with different numbers of bins. As 
measurements only exist for 8 pairs of bins in the case of .Kg7r + 7T~ and for 2, 3 and 
4 pairs of bins in the case of K§K + K~ , the study instead uses input values of (q, 



Si) and Tj as calculated from the model in Ref. 25 , and these are assumed to be 
known perfectly. The free parameters of the fit are xp, Vd, t cp-, a cp, T and the total 
number of decays. For K^7t + it~ experiments are performed with 8, 16, 32 and 50 
pairs of equal-phase bins, and for K§K + K~ studies are made with 4, 8, 16 and 32 
pairs of bins. 

The precision determined for each scenario is plotted in Figs. [2] and [3] for 100M 
K^tt + tt~ and 20M K§K + K~ decays. In general the precision improves with the 
number of bins. In order to quantify this trend the following function has been fit to 
each variable v: 

a(v) = a+^ (16) 

where N is the number of Dalitz plot bin pairs. The asymptotic parameter a gives 
the statistical precision that is achievable in the unbinned scenario. The values of a 
and b for each mixing and CP violation parameter are shown in Table |6j 

Although improved precision is in principle achievable by increasing the number of 
bins from the values for which measurements are currently available, this improvement 
is in all cases rather modest, even when going to an infinite number of bins (~ 10% 
for most parameters). This indicates that the model-independent analysis with 8 
pairs of bins for K^7t + tt~ , and 4 pairs of bins for K$K + K~ , has a statistical precision 
which is only marginally inferior to that of the model-dependent, unbinned analysis. 
Furthermore there can be no significant improvement in sensitivity by changing the 
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Figure 2: Dependence of the fit uncertainty for the mixing and CP violation pa- 
rameters against number of Dalitz bin pairs for 100M K^n + ir~ decays. The input 
parameters and fit conditions are described in the text. Superimposed is a fit to 
Eqn. 16 (solid line) and the asymptotic value of the fit (dotted). 
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Number of bin pairs Number of bin pairs 




Number of bin pairs Number of bin pairs 



Figure 3: Dependence of the fit uncertainty for the mixing and CP violation pa- 
rameters against number of Dalitz bin pairs for 20M K®K + K~ decays. The input 
parameters and fit conditions are described in the text. Superimposed is a fit to 
Eqn. 16 (solid line) and the asymptotic value of the fit (dotted). 
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Table 6: Results for the parameters a and b, defined in Eqn. 16, that describe the 
dependence of the precision of the mixing and CP parameters on the number of Dalitz 
plot bins. 







+ 7T- 






Parameter 


a 


b 


a 


b 


a{x D ){%) 

&(r C p) 
<t(«cp)(°) 


0.0130(2) 
0.0132(2) 
0.0121(2) 
0.754(9) 


0.0122(25) 
0.0063(25) 
0.0059(25) 
0.668(130) 


0.0373(3) 
0.0218(2) 
0.0283(2) 
1.789(14) 


0.0350(21) 
0.0060(11) 
0.0014(15) 
2.091(115) 



binning scheme in the analysis from the equal-interval strong-phase binning which 
has been assumed for the studies presented here. In practice the gain in precision 
when increasing the number of bins will presumably be even more modest than Figs. [2] 
and [3] suggest, as the systematic uncertainties associated with the knowledge of (q, 
Si), not included in this study, will grow as the bin sizes become smaller. 

5 Discussion and summary 

As already demonstrated by existing measurements, a time-dependent analysis of 
K^h + h~ decays is a powerful method for studying mixing and to probe for indirect 
CP violation in the D°-D° system. Observation of CP violation in charm mixing 
related observables at a level beyond the very small value predicted in the Standard 
Model would be a strong indication of New Physics contributions. A binned, model- 
independent fit is an interesting alternative analysis strategy to the unbinned model- 
dependent approach which has so far been pursued. An extensive programme of 
simulation studies has been performed to assess the strengths and limitations of the 
binned analysis. From these studies three main conclusions emerge: 

1. A binned analysis is a very attractive approach for extracting mixing and CP 
violation information from K®h + h~ decays. By making use of the measurements 
of the (cj, si) parameters available from CLEO-c it is possible to obtain a 
statistical precision on xp, yo, rep and acp that is only marginally inferior to 
that available from an unbinned analysis. Most importantly, the binned analysis 
has no systematic error arising from the use of an amplitude model, in contrast 
to the unbinned approach where this source of uncertainty is unavoidable. 

2. The existing CLEO-c measurements of (q, Sj), performed with only 0.8 fb _1 of 
■0(3770) data are adequate for analysing the size of sample that can be expected 
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at the current LHCb experiment. Looking further forward to the extremely large 
samples that may be accumulated at the LHCb upgrade or Super- B or Belle-II 
it is very unlikely that the uncertainties on c, and Sj will ever become limiting 
in the measurement of rep and acp, which are the principal parameters of 
interest. Even with 100M K^tt + ti^ decays the expected statistical uncertainty 
on the CP parameters is only around half that of the systematic arising from 
the CLEO-c measurement errors on q and Sj. BES-III already has a ^(3770) 
sample that is four times larger than that analysed by CLEO-c, and this is set 
to grow. Analysis of these data will allow for the (q, si) systematic uncertainty 
to become sub-dominant, even for enormous samples of K^h + h~ decays. Plans 
exist for future projects, such as a Super-5 threshold run, or a Novorsibirsk 
r— charm factory, that will provide a ^(3770) data set one-to-two orders of 
magnitude larger than that accumulated at BES-III. Although welcome for 
many reasons, these initiatives are not essential for improving the knowledge of 
the CP violation parameters from the analysis of flavour-tagged K§h + h~ data. 
On the other hand, the (c$, Sj) uncertainties will more rapidly become limiting 
for the mixing parameters xp> and yp, and here larger quantum-correlated D D 
samples will be useful. 

The bin granularity and partitioning schemes for which the CLEO-c measure- 
ments have been performed look well suited for D°-D° mixing and CP violation 
studies with all current and future data sets. 

3. With very large K§h + h~ samples an alternative fit strategy becomes possible 
which will benefit the measurement of the unitarity triangle angle 7/^3 using 
B~ — > D(Ksh + h~)K~ decays. The induced systematic uncertainty on 7/^3 
arising from the CLEO-c measurements of (q, Si) is estimated to be 2.1° for the 



equal phase binning 18 . This uncertainty is small compared to the anticipated 
statistical precision of the current LHCb experiment, but may become limiting 
at the upgrade, where a measurement error of ~ 1° is expected. A similar 
precision is foreseen at the Super-5 experiment. It is therefore desirable to 
reduce the uncertainties on (q, Si) by a factor of two or better with respect to 
the CLEO-c analysis. 

Such an improvement may be possible from analysis of the BES-III data alone. 
However, a complementary measurement of similar precision can be performed 
using flavour-tagged K$h + h~ data. In this approach external constraints are 
placed on i D and yp, and the fit is used to determine (q, Sj), as well as r C p and 
acp- Assuming 100M _Kg7r + 7r~ events, and a knowledge of xp, and yp> which is a 
factor four times better than at present, then the uncertainties on q and can 
be decreased to around 60% (on average) of their current values. This strategy 
therefore appears promising, although it requires improved knowledge of xp, 
and yp, which must come from independent data sets and/or mixing analyses 
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based on other decay modes. 

In conclusion, the model-independent binned fit has many attractive features to 
recommend it as an analysis strategy for K®h + h~ decays, both with the data sets 
currently available and the very large samples that will be collected over the coming 
decade. First results using this approach are eagerly awaited. 

Acknowledgments 

We acknowledge useful discussions with Anton Poluektov and Jim Libby. We are 
grateful for support from the STFC, United Kingdom. 

References 

[1] M. Staric et al. (Belle Collaboration), Phys. Rev. Lett. 98 (2007) 211803; 
L.M. Zhang et al. (Belle Collaboration), Phys. Rev. Lett. 96 (2006) 151801. 

[2] B. Aubert et al. (BaBar Collaboration), Phys. Rev. D 80 (2009) 071103 
B. Aubert et al. (BaBar Collaboration), Phys. Rev. Lett. 103 (2009) 211801 
B. Aubert et al. (BaBar Collaboration), Phys. Rev. D 78 (2008) 011105(R) 
B. Aubert et al. (BaBar Collaboration), Phys. Rev. Lett. 98 (2007) 211802. 

[3] T. Aaltonen et al. (CDF Collaboration), Phys. Rev. Lett. 100 (2008) 121802. 

[4] L.M. Zhang et al. (Belle Collaboration), Phys. Rev. Lett. 99 (2007) 131803. 

[5] P. del Amo Sanchez et al. (BaBar Collaboration), Phys. Rev. Lett. 105 (2010) 
081803. 

[6] A. Zupanc et al. (Belle Collaboration), Phys. Rev. D 80 (2009) 052006. 
[7] R. Aaij et al. (LHCb Collaboration), JHEP 04 (2012) 129. 
[8] R. Aaij et al. (LHCb Collaboration), Phys. Rev. Lett. 108 (2012) 111602. 
[9] CDF Collaboration, |arXiv: 1207.2158] [hep-ex]. 



[10] D. Asner et al. (Heavy Flavor Averaging Group), arXiv: 1010. 1589 [hep-ex]; online 



updates at http://www.slac.stanford.edu/xorg/hfag ; quoted numerical results 
are from December 2011. 

[11] S. Bianco, F.L. Fabbri, D. Benson, I. Bigi, Riv. Nuovov. Cim. 26N7 (2003) 1. 

[12] D.M. Asner et al. (CLEO Collaboration), Phys. Rev. D 72 (2005) 012001. 



17 



T. Aushev et al, |arXiv:1002.5012| [hep-ex]. 



B. Meadows et al, |arXiv: 11 09.5028 [hep-ex]. 
LHCb Collaboration, CERN-LHCC-2011-001. 
LHCb Collaboration, CERN-LHCC-2012-007. 

A. Bondar, A. Poluektov and V. Vorobiev, Phys. Rev. D 82 (2010) 034033. 

J. Libby et al. (CLEO Collaboration), Phys. Rev. D 82 (2010) 112006. 

S. Malde and G. Wilkinson, Phys. Lett. B 701 (2011) 353. 

J.L. Rosner et al. (CLEO Collaboration), Phys. Rev. Lett. 100 (2008) 221801; 
D.M. Asner et al. (CLEO Collaboration), Phys. Rev. D 78 (2008) 012001; N. 
Lowrey et al. (CLEO Collaboration), Phys. Rev. D 80, (2009) 031105(R); J. 
Insler et al. (CLEO Collaboration), Phys. Rev. D 85 (2012) 092016. 

A. Giri, Y. Grossman, A. Soffer and J. Zupan, Phys. Rev. D 68 (2003) 054018. 

A. Bondar and A. Poluektov, Eur. Phys. J. C 47 (2006) 347; A. Bondar and 

A. Poluektov, Eur. Phys. J. C 55 (2008) 51. 

B. Aubert et al. (BaBar Collaboration), Phys. Rev. D 78 (2009) 034023. 

A. Poluektov et al. (Belle Collaboration), Phys. Rev. D 81 (2010) 112002. 

P. del Amo Sanchez et al. (BaBar Collaboration), Phys. Rev. Lett. 105 (2010) 
121801. 

L. Dong, 'Recent Results from BESIII', EPS HEP2011, Grenoble, France, July 
2011. 

A. Bondar and E. Levichev, 'Physics case of tau-charm factory and Novosibirsk 
plan', ECFA meeting, CERN, November 2010. 



18 



